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ABSTRACT 

We use a preliminary version of our "NextGen" grid of cool star model atmo- 
spheres to compute synthetic line profiles which fit high resolution Keck spectra of the 
cool M dwarf VB10 satisfactorily well. We show that the parameters derived from the 
Keck data are consistent with the parameters derived from lower resolution spectra 
with larger wavelength coverage. We discuss the treatment of van der Waals broad- 
ening in cool, molecular (mostly H 2 ) dominated stellar atmospheres. The line profiles 
are dominated by van der Waals pressure broadening and are a sensitive indicator 
for the gravity and metallicity. Therefore, the high-resolution Keck spectra are use- 
ful for determining the parameters of M dwarfs. There is some ambiguity between 
the metallicity and gravity. For VB10, we find from the high-resolution spectra that 
5.0 < log(g) < 5.5 and < [^-] < +0.5 for an adopted fixed effective temperature of 
2700 K (Schweitzer, 1995), which is consistent with recent interior calculations (e.g. 
Baraffe et al., 1995). 
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1 INTRODUCTION 

M dwarfs are among the faintest and coolest stellar objects. 
Their spectra are dominated by molecular band absorption. 
In the optical region the major opacity sources are the TiO 
bands, which produce a pseudo continuum with the atomic 
lines superposed. These atomic lines show very broad wings 
due to van der Waals broadening, whereas the molecular 
lines show only relatively small damping wings. We have 
recently obtained a high resolution optical spectrum of the 
dM8e star VB10 (Gl 752B) with the HIRES echelle on the 
Keck telescope. The collection and reduction of the data are 
as in Basri & Marcy (1995). In this paper, we use this spec- 
trum to assert our treatment of the van der Waals line broad- 
ening in M dwarf atmospheres and to estimate the gravity 
and the chemical composition of VB10. VB10 is known to 
be a chromospherically active flare star. However, in this pa- 



Based on observations obtained at the W.M. Keck Observatory, 
which is operated jointly by the University of California and the 
California Institute of Technology 



per we will neglect the chromosphere and do not attempt to 
model chromospheric spectral features, such as core emission 
observed in resonance lines. 

Our calculations are based on a preliminary version of 
our latest grid of model atmospheres (the "NextGen" or 
"version 5" grid) for cool dwarf stars and version 6.2 of the 
generalized stellar atmosphere code phoenix. An important 
improvement over the previous generation of this grid (Al- 
lard & Hauschildt, 1995, hereafter, AH95) are larger and 
more reliable molecular line lists which allow us to com- 
pute detailed molecular line profiles. A complete descrip- 
tion of all changes and improvements over the AH95 version 
of the models will be given in a subsequent paper (Allard, 
Hauschildt, & Schweitzer, 1995, hereafter, AHS96). Note 
that we use metallicities scaled from the solar values of An- 
ders & Grevesse (1989) and we use the notation 



log 



/H 



to specify abundances 
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A previous analysis based on low resolution spectra by 
Brett (1995) yielded only rough values for T cti =2400-2600 
K, log (3) w 5.0 and [^1 « 0.0. Interior calculations (e.g. 
Burrows et al., 1993) require for solar M dwarfs in the 
effective temperture range of « 2700 K gravities of log(<?) « 
5.3. 

In the next section we will describe our theoretical ap- 
proach to the problem. We will discuss the basic approxima- 
tion we use as well as their numerical realization. In section 
3 we present the results of our modeling. We discuss the sen- 
sitivities of the profiles and try to establish error limits. In 
section 4 we compare the observed high resolution spectra 
of VB10 to our model spectra. 



the equation for Ce which leads to 
C° = 1.01 x 10~ 32 (Z + l) 2 



El 



(E - Ei) 2 (E- 



r 6 -li 
cm s 



(3) 



where a p is the polarizability of the perturber, Z the charge 
of the absorber, E the ionization potential, Ei the lower 
and E u the upper level excitation energy of the absorber, 
qh the polarizabilty of Hydrogen, and Eh = 13.6eV. We 
include the most abundant perturbers in M dwarfs and their 
polarizabilities as given in Weast (1988) and listed in Tab. 1 
to calculate the total Van der Waals damping constant: 
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2 THE TREATMENT OF LINE BROADENING 
2.1 The line profiles 

We include all common line broadening mechanisms, i.e., 
thermal broadening, micro-turbulence, natural broadening 
and pressure broadening, both in the model calculations 
and synthetic spectra. The resulting Voigt profiles V(u,a), 
where u = AA/AAd and a = (A 2 7/47tc)/AA.d (Ad is the 
Doppler width, 7 the Lorentz width, AA the distance from 
the line center and A the wavelength), are calculated as a 
convolution of a Gauss and a Lorentz profile using stan- 
dard approximations (e.g. Baschek & Scholz, 1982). The gas 
temperatures in M dwarf atmospheres are not high enough 
to sustain a significant amount of ionization in the atmo- 
sphere. The electron and proton densities are, therefore, 
much smaller than the densities of the most important neu- 
tral and molecular species. Consequently, the contribution 
of Stark broadening to the total damping constant is very 
small, even in stars with very low metallicities. We include a 
microturbulence of £ = 2 km sec -1 in addition to the ther- 
mal speed of the atoms and molecules. However, the to- 
tal thermal plus microturbulent line widths are for most 
lines much smaller than the line width due to van der Waals 
(vdW) broadening. In the following paragraphs we will de- 
scribe in some detail the way we estimate the vdW damp- 
ing constants and how their choice influence the computed 
profiles of both atomic and molecular lines. The interaction 
between two different, unpolarized, neutral particles is de- 
scribed by the van der Waals interaction 



Av = 



C 6 



(1) 



where Av is the resulting frequency shift due to one interac- 
tion, r the relative distance between the interacting particles 
and Ce the interaction constant for the vdW interaction. 
Within the impact or static approximation vdW damping 
leads to a Lorentz profile with a full width half maximum 
damping constant 7vdw (cgs-units) 



7vdw = 17 • C, 



2/ \ 3 / 5 N v . 



(2) 



Here v is the relative velocity between perturber and ab- 
sorber and N p the number density of the perturber. 

We base the calculation of Ce on Unsold's hydrogenic 
approximation (Unsold, 1955). However, we explicitly ac- 
count for the different polarizabilities of each perturber in 



where we sum over all perturbers considered. Due to the 
lack of availability of any vdW broadening mechanismns for 
molecules in this temperature range, we also use this ap- 
proximation for the molecular lines. 

Earlier investigations (Weidemann, 1955) showed that 
the values as calculated by Eq. (3) are in good agreement 
with observed line widths for alkali metals but not for other 
elements, such as iron (Kusch, 1958). This has lead to the 
introduction of correction factors to the 'classical' formula. 
Therefore, we include a correction factor for C6 of 



C 6 = CT 1 x C2 = 10 1 ' 8 x C°e 



(5) 



as described in Wehrse & Liebert (1980) for non-alkali like 
species. For alkali-metals and ions with alkali-like electron 
structure we omit any correction factor. This will be dis- 
cussed and investigated in greater detail in sections 4 and 5. 
However, as these results will show, the use of this correction 
factor remains unclear for the non-alkali like species. 



2.2 The numerical treatment of vdW broadening 

In order to calculate the emerging spectrum we have to cal- 
culate a Voigt profile at every wavelength point for every 
line at every depth point. To limit the computation time we 
make the following simplification. 

Before calculating the absorption coefficient, a line se- 
lection procedure decides which lines are calculated with 
Voigt profiles. We follow the line selection procedure de- 
scribed by AH95. Only lines that are stronger than a cer- 
tain threshold relative to the continuum are treated with 
Voigt profiles. This is decided by comparing the absorption 
coefficient at the line center Ki with the corresponding con- 
tinuous absorption coefficient k c . Only if the ratio V = ki / k c 
for at least one of three representative depth points is larger 
than a preset threshold value, a line is treated with a Voigt 
profile. Otherwise the line is considered so weak that to a 
good approximation it can be treated as a Gauss profile in 
order to save computation time. The number of considered 
atomic lines is much smaller than the number of molecular 
lines. This allows us to set T for atomic lines very low in 
order to treat as many lines as possible with Voigt profiles. 
The number of molecular lines is about two orders of mag- 
nitude larger but since the atmosphere structure does not 
change significantly (as test calculations showed) with the 
inclusion of more Voigt profiles for molecular lines we set the 
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Table 1. The perturbers included in the calculations and their polarizabilities. The percentages are taken from models with log(g) = 5.0, 
[Ml = ^ anc ' t ' le res P ec ti ve T e ff. We chose the optical depth r stt j = 1 (dr st( j = nds where k is the absorption coefficient at a standard 
wavelength of 1.2^m (AH95) and ds the geometric depth) as an example for an average line forming depth. The percentages are rounded 
on the last figure. 



Perturber 


a p in l(r 24 cm 3 


Fraction of P g 
T eS =2000 


;as in per cent at 
' T eft =2700 


fstd = 1 for 
T cH =3500 
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0.666793 


1.2 


14.9 


66.6 


He 


0.204956 


16.2 


15.1 


10.7 


Ne 


0.3956 ±0.1% 


0.02 


0.02 


0.01 


Fc 


8.4 ± 25% 


0.01 


0.01 f 


s 0.005 


H 2 


0.806 ± 0.5% 


82.4 


69.8 


22.4 


CO 


1.95 ±0.5% 


0.06 


0.06 


0.04 


H a O 


1.45 ±0.5% 


0.07 


0.06 


0.01 


N 2 


1.7403 ± 0.5% 


0.01 


0.01 f 


s 0.004 



Figure 1. The calculated Nal AA8183.8195 doublet in vac- 
uum with different broadening mechanisms (log(g) = 5.0, 
T c ff = 2700.K", = 0.0). From top to bottom : Gauss pro- 

file, only H and He as perturbers, H, He and H 2 as perturbers, 
perturbers of Tab. 1 without the correction factor (this plot is 
lying exactly on the previous one), perturbers of Tab. 1 and the 
correction factor of 10 18 artificially applied. 

threshold for molecular lines for iterations to V — 100. For 
high resolution spectra we use lower values (usually r w 1) 
although we found that the spectra do not change signifi- 
cantly by increasing the number of lines with Voigt profiles. 

The spectra presented throughout this paper used for 
molecular lines F = 5 when we examined atomic lines and 
r = 0.1 when we examined molecular lines. All the spectra 
are furthermore based on self-consistent model strucutres 
with thier respective parameters. 



3 RESULTS 

3.1 Synthetic line profiles 

To demonstrate the various effects influencing the line pro- 
files we will use the Nal AA8183, 8195 subordinate doublet 
and the Kl AA7665, 7699 resonance doublet, since they are 
available in the Keck-spectra (see Sec. 4). Note that the syn- 
thetic line profiles in this section are shown at their vacuum 
wavelengths. 

Eq. (2) shows that the vdW damping constant is a lin- 
ear function of the perturber density. This implies that the 
major contribution to the total vdW broadening will be due 
to the effect by the most abundant perturbers, namely Hi, 
He I and H 2 , because the polarizabilities change much less 
than the relative concentrations of the perturbing species 
(cf. Tab. 1). In Fig. 1 we show the influence of H 2 on the 
profile of the subordinate Nal AA8183, 8195 doublet. All the 
other perturbers listed in Tab. 1 do not contribute signifi- 
cantly to the line width, the profile is identical with the one 
which is only H-He-H 2 broadened. This is valid even for 
metal abundances higher than solar. 

In addition, Fig. 1 shows the effect of the correction 
factor introduced in Eq. (5) by artificially applying it to the 
alkali metal sodium. In this case the correction factor was 
only applied to demonstrate the uncertainty which still is 



Figure 4. The calculated Nal AA8183, 8195 doublet in vacuum 
with different broadening methods (log(g) = 5.0, T c ff = 2700-K, 
[Ml = o.O). From top to bottom : our method without the cor- 
rection factor, using Kurucz's damping constant per Hydrogen 
atom at 10 000 K, our method and the correction factor of 10 18 
artificially applied, the old method. See Sec. 3.2 for details. 

present in the vdW damping constant. A final decision on 
the validity of such a factor has to be made by comparing 
the synthetic profiles with observations. 

Changes of the absolute number density of each per- 
turber will change the damping constant. We consider two 
ways to change the perturber densities: changing the grav- 
ity or changing the metallicity. Higher surface gravities will 
increase the total pressure in every layer of the atmosphere. 
This will also increase the number density of every perturber 
for constant metallicity. An example of the gravity effect on 
the line widths is shown in Fig. 2 which illustrates the influ- 
ence of log(g) on the profile of the subordinate Nal doublet 
as well on the K I AA7665, 7699 resonance doublet for model 
atmospheres with T c ff = 2700 K, solar metallicity and a dif- 
ference in the gravities of one order of magnitude. 

On the other hand, an increase in metallicity at con- 
stant gravity will decrease both the total pressure and the 
number density of the main perturbers, Hi, He I and H 2 . 
It furthermore will bind hydrogen in hydrides and water 
vapour. Therefore, all these effects will make the line wings 
narrower. We demonstrate this in Fig. 3, with the same two 
doublets as above. The only difference in the models shown 
is half an order of magnitude in the metallicity. 

This introduces the problem that small changes in 
metallicity can be compensated for by changes in the grav- 
ity, thus making the determination of both parameters from 
the wings of strong lines alone questionable. This means 
that abundance analyses of cool dwarfs in general must be 
performed with caution and simple techniques may lead to 
errors. This will be investigated in much more detail in sub- 
sequent work. 

3.2 Comparison to other methods 

To demonstrate the improvement brought by the present 
handling, we will briefly review other commonly used meth- 
ods: (1) the Unsold formula without perturber correction 
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Figure 2. The calculated Nai AA8183, 8195 and Ki AA7665, 
(lower line)(T cff = , [M] = 0.0) 



Figure 3. The calculated Nai AA8183, 8195 and Ki AA7665, 7699 
line) (T eS = 2700E-, log(g) = 5.5) 

used by Allard (1990) and AH95, and (2) using the precal- 
culated damping constants per hydrogen atom at 10 000 K 
provided on CD-ROM No.l by Kurucz (1994). 

The method (1) used an average broadening mechanism 
by taking the total particle density instead of the density of 
every perturber and by not distinguishing between the dif- 
ferent perturbers, i. e., not taking into account the different 
polarizabilities and reduced masses, latter influencing the 
relative velocity in Eq. (2). In addition, it used the correc- 
tion factor for all lines. 

In Fig. 4 we compared now the same Nai AA8183, 8195 
profiles as above calculated with the old method with one 
calculated with the new method but the correction factor 
artificially applied for comparison. As can be seen there are 
clear differences between the two methods. The new method 
results in narrower lines since it accounts for the different 
polarizabilities and reduced masses which decrease Ca in 
Eq. (3) or 7 v< jw in Eq. (2) compared to a C\ and 7 v< jw 
calculated as described above and, therefore, also decreases 
the damping constants of Eq. (4). 

In Fig 4 we also plot the profiles using the damping 
constants provided by Kurucz. We find only small differ- 
ences between Kurucz's vdW damping constant and ours, 
if we omit the correction factor. In test calculations we also 
compared the use of natural damping constants provided by 
Kurucz with those calculated with the classical approxima- 
tion. Since the total line width in M dwarfs is dominated by 
vdW broadening, it makes little observable difference which 
treatment of natural damping is employed to model their 
spectra. We adopt our own methods in the following discus- 
sion. 

We furthermore compared directly the damping con- 
stants obtained with our approximation with the corre- 
sponding values from Kurucz's list and found only small 
differences for the vast majority of the lines above 6000 A. 
The largest differences arise from the correction factor we 
included for non-alkali species. For the UV region the differ- 
ences are larger. Here the hydrogenic approximation for the 
absorber we used differs substantially from the method Ku- 
rucz used. This does not affect the spectra we present here. 
In general, the photospheric UV flux from M dwarfs is negli- 
gible. However, a more thorough investigation must be per- 
formed for stars that are significantly hotter than M stars. 

3.3 Uncertainties in the broadening mechanism of 
molecular lines 

As can be seen from Figs. 5 and 11 the molecular lines are 
much narrower than the atomic lines, both in the observa- 
tions (cf. Sec. 4.6) as well as in our models. This is due to 
the fact that each molecular line is very weak compared to 
the strong atomic lines. Therefore only the Gauss cores re- 
main visible in the spectra and the line width is only the line 



doublets in vacuum for log(<;) = 4.5 (upper line) and log(<;) = 5.5 



in vacuum for |^ jj- J = 0.5 (upper line) and = 0.0(lower 



width of the Gauss core. Nevertheless, the Voigt profiles need 
to be included since they lift the Gauss cores significantly 
without changing the line width as corresponding compar- 
isons showed. In addition, the molecular lines lie so densely 
together that only their cores can remain visible; any wing 
contribution will effectively be blocked by the core of neigh- 
boring molecular lines. This means that the total width of 
molecular features is not as much dominated by the wings 
of the vdW broadening as for the atomic lines. 

In order to estimate the error resulting from our approx- 
imations, we performed test calculations in which we omit- 
ted the correction factor from Eq. (5) for molecules which 
results in a change in the interaction constant by nearly 2 
orders of magnitude. In contrast to the atomic lines (cf. Sec. 
3.1 and Fig. 1), we found only small differences in the over- 
all appearance of molecular lines in the optical (see Fig. 5), 
since only the line wings change and they are completely 
blended in the spectrum. 

The most important result, however, was the influence 
of rotational broadening on the appearance of the molecular 
lines which can be seen in Fig. 5 as well. As we applied rota- 
tional broadening on our test spectra we found that even a 
small rotational velocity will broaden the lines so much that 
the uncertainties in the vdW broadening become unimpor- 
tant. Only the cores of the stronger molecular lines remain 
visible. Also the influence of the gravity and metallicity be- 
come less important. We furthermore found that the effect of 
rotational broadening rises continuously with increasing ro- 
tational velocity. As a consequence, the molecular line 'haze' 
reveals a possibility to accurately measure the rotational ve- 
locity as it is very sensitive to it (see Sec. 4.6). 



4 COMPARISON OF VB10 

HIGH-RESOLUTION KECK SPECTRA 
WITH SYNTHETIC SPECTRA 

Our observation of VB10 was obtained with the HIRES 
echelle on the Keck telescope during a run on 1995 March 
12 under clear conditions. Our exposure time was 30 min- 
utes with 0.8 arcsec seeing. The instrumental setup and data 
reduction were very similar to those described by Basri & 
Marcy (1995). The wavelength setting was slightly modi- 
fied to include the subordinate Nai AA8183,8195 doublet, 
the Ki AA7665, 7699 resonance doublet, the Cal A6573 res- 
onance line, the Bal A7911 resonance line and the Fe I A7913, 
Fcl A8662 and Fel A8689 subordinate lines. Inclusion of 
Nai meant that we could not also observe the Rbl line 
discussed by Basri & Marcy (1995). 

We consider models with any combination of \og(g) = 
4.0, 4.5, 5.0, 5.5 and scaled solar metallicities = 
—0.5, 0.0, +0.5. For the purpose of this paper, we adopt a 
fixed effective temperature of 2700 K for VB10, derived from 
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Figure 5. Calculated molecular lines for log(g) = 5.0, I jj- J = +0.5, T e g = 2600i^. Dotted : Without the artificial correction factor 
of 10 1,8 , Solid : With the correction factor applied. Left panel : Not rotationally broadened. Right panel : rotationally broadened with 
lOkms" 1 . 



fitting low-resolution spectra obtained by Kirkpatrick et al. 
(1993) and Jones et al. (1994) to low-resolution models (see 
Schweitzer, 1995). It is not the purpose of this paper to in- 
vestigate the influence of the effective temperature on the 
line profile. This will be done in subsequent work where the 
parameters will be verified or revised, if necessary. A rota- 
tional velocity of 8kms -1 was adopted, using the method 
of Basri & Marcy (1995). We reconsider this in section 4.6. 
Due to the complexity of the spectra we performed the fits 
and decided on their quality by eye. 

4.1 The sodium doublet 

The Nal AA8183,8195 doublet lines are subordinate tran- 
sitions and, therefore, form mainly in the photosphere and 
should not be significantly affected by the chromosphere of 
VB10. Our best fits are the two models with log(g) = 5.0 and 
[f ] = 0.0 and with \og(g) = 5.5 and [f ] = +0.5 shown in 
Fig. 6. There are no significant differences in the appear- 
ance of the two model spectra, although the parameters 
differ substantially. This demonstrates the effect mentioned 
above, namely that an increase in metallicity decreases the 
number density of the major perturbers, Hi, He I and H2, 
and thus partly cancels the effect of increasing the grav- 
ity, in particular for strong lines. This shows that one has 
to be very careful when analyzing M dwarf spectra in de- 
tail. In our case here, we found good agreement between 
observation and models for parameters within the intervals 
5.0 < log(g) < 5.5 and 0.0 < [f ] < +0.5. We therefore use 
this as an error limit on our analyses. For any parameter 
combination beyond these intervals we find no reasonable 
agreement. 

The synthetic line profiles show a deeper core than the 
observed lines. We calculated several models with a different 
microturbulence in order to "fill in" the core but we found no 
significant improvement to the fit. The remaining difference 
is most likely due to NLTE effects which might make line 
cores less deep. This will be investigated in detail in future 
work (see e.g., Hauschildt et al. (1996)). NLTE effects will 
have to be considered as they are highly non-linear and are 
very likely to affect the cores of strong lines. 

We also like to point out the numerous weak lines in 
the wings of the atomic lines. These are mostly TiO lines 
which are also visible in the observed spectrum and are not 
to be considered as noise (see also below, Sec. 4.6). 

4.2 The potassium doublet 

The Ki AA7665, 7699 lines show a clear core reversal (see 
Fig. 7). This is due to the chromospheric activity of VB10. 
The doublet is a resonance transition and thus will be 
strongly influenced by chromospheric activity. Since we are 
not treating the chromospheric effects here, our models can 
only reproduce the line wings which are formed deeper in 
the photosphere. 

As in the case of the sodium lines, both models fit 



equally well and a distinction between the two parameter 
sets of log(g) = 5.0 and [^1 = 0.0 and of log(g) = 5.5 and 
[M] = + o.5 is not readily possible (see Fig. 7). The fit it- 
self is not as good as in the case of the sodium doublet as 
the models show many more molecular lines. This spectral 
region is partially contaminated by the telluric O2 A bands 
present between 7595A and 7680A which we did not remove 
from the observed data. They are affecting the observations 
by decreasing the continuum flux to a not readily known 
extend. In addition, this part of the spectrum is dominated 
by the e-band of TiO for which the line data are of lower 
quality (AHS96). Nevertheless, the lines are so broad, that 
we still can fit them as shown in Fig. 7. 

4.3 The calcium line 

The Ca I A6573 resonance line shows also chromospheric fea- 
tures as can be seen in Fig. 8. The line is weak and very 
much filled up and the line is hardly detectable without a 
specific identification. Calcium, in contrast to the two ele- 
ments discussed above, is not an alkali metal. This means 
that the vdW damping constant is calculated including the 
correction factor introduced in Eq. (5). Although it is hard 
to compare this weak line with models, we find that the cor- 
rection factor is necessary in order to reproduce the observed 
line width with the same model parameters as used for the 
sodium and potassium lines. Without C6° rr the theoretical 
line profile of the Cal lines is not broad enough, however, 
we cannot say how accurate our value for Ce° rr is. In order 
to make a more accurate determination of this factor, wc 
will have to investigate much stronger non-alkali lines with 
clearer and broader wings. This will be done in a subsequent 
analysis. 

This line cannot be used to determine directly the pa- 
rameters of VB10. But it can be used to confirm them. As 
in the case of Sodium and Potassium above, the two mod- 
els with log(gr) = 5.0 and [MJ — fj.O and log(g) = 5.5 and 
[If] ~ +0.5 produce equally good results. 

4.4 The barium line and the iron line at 7913 A 

The fit to the Bal A7911 resonance line and the Fel A7913 
line (which has 0.86 eV excitation energy) is shown in Fig. 
9. Both lines are very weak but still visible and we find very 
good agreement between observations and our two models 
with log(ff) = 5.0 and [M] = 0.0 and log(ff) = 5.5 and 
[M] = +0.5. As in the case of calcium above, we only used 
these lines to confirm our range of fitting parameters, since 
the line widths are not dominated by the wings of the pro- 
files. They are dominated by their Gauss cores but the Voigt 
profiles need to be included in order to lift the cores. 

Both metals, barium and iron, are non-alakli metals 
as calcium above and, therefore, treated with Cg orr from 
Eq. (5). Although the lines are narrow, we find it also here 
necessary to include such a correction factor in order to fit 
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Figure 6. The fit to the Nal AA8183, 8195 doublet. Solid : Observed. Dotted : T cB = 2700^, log(g) = 5.0, [ M ] = 0.0 (left) and 
T cS = 2700K", log(g) = 5.5, [Mj = +0.5 (right). Both models with Skms" 1 rotationally broadened. 

Figure 7. The fit to the Kl AA7665, 7699 doublet. Solid : Observed. Dotted : T off = 2700/f , log(g) = 5.0, [M] = 0.0 (left) and 
T cff = 2700A", log(g) = 5.5, [— ] = +0.5 (right). Both models with 8kms _1 rotationally broadened. 



all lines consistently. However, we cannot make any accu- 
rate determinations of Cf° IT , which still has to be done with 
stronger lines. 

We also want to point out the very good fit of the sur- 
rounding molecular lines (see also Sec. (4.6)). 

4.5 The iron lines at 8662 and 8689 A 

The two subordinate Fe I A8662 and Fe I A8689 lines visible 
in the observations are shown in Fig. 10. The lines are not 
very strong, yet they can be clearly identified. This spectral 
region is strongly contaminated by terrestial absorption fea- 
tures which reduce the total flux and suppress the molecular 
pseudo continuum. Therefore, we did not try to fit the lines 
in detail. Instead we show them with an offset to demon- 
strate the resemblance of the line shapes of the two Fel 
lines. 

We cannot use these lines to determine the parameters 
of VB10. As in the cases above we used models with \og(g) = 
5.0 and [f] = 0.0 and log(ff) = 5.5 and [f] = +0.5 and 
found consistency within the given constraints. Again, cal- 
culating the vdW damping constant including Cg orr (see also 
Sec. (4.4), results in consistent synthetic spectra. But due 
to the uncertainties in the observed spectrum we did not try 
to determine the accuracy of Cg orr here either. 

We want to point out the 'edge' at approximately 8671 
A in our models, which is due to an abrupt end of a VO band. 
In our models we can treat VO only with the JOLA method 
(AH95, AHS96) since no line list is yet available. Since the 
JOLA method will produce these edges, it is not suitable 
for computing high-resolution spectra. We point out that 
the Fe I A8662 line is slightly contaminated by a Ca II A8662 
line. 

4.6 Molecular lines 

As mentioned before, the high-resolution spectra are full 
with a large number of weak molecular lines. For TiO, we use 
the semi-empirical line list of J0rgensen (1994), which will 
not reproduce most of the TiO lines at their exact wave- 
length or their proper gf-value. But nevertheless many of 
the lines agree very well with observations. We found an ex- 
cellent agreement with models of the two parameter sets 
of log(ff) = 5.0 and [f] = 0.0 and log(g) = 5.5 and 
[Ml = +0.5 as can be seen in Fig. 11. Because of the 
strong influence of rotational broadening, we only use the 
molecular lines to confirm the parameters and not to derive 
them directly. On the other hand, we can derive the rota- 
tional velocity a second time by using its strong influence 
on the profile of the molecular lines. We found that rota- 
tional velocities between 7 and 9kms _1 produce very simi- 
lar fits, whereas applying values beyond this interval result 



in worse fits. Therefore, we can confirm a rotational velocity 
of 8±lkms-\ 



5 SUMMARY AND CONCLUSIONS 

We found that the line widths caused by pressure broadening 
due to van der Waals damping are very sensitive to both 
the gravity and the metallicity. Therefore fitting the lines 
of high resolution spectra provides a very important tool 
to determine these two parameters. Nevertheless, since both 
parameters influence the line widths in the similar ways, one 
has to be very careful with the analysis. 

From the fits to the lines described above we estimate 
that VB10 has a gravity in the range 5.0 < log(g) < 5.5 
and a metallicity 0.0 < [— ] < 0.5 based on line broadening 
considerations alone with an adopted effective temperature 
of T cff = 2700K (Schweitzer, 1995). We will investigate the 
influence of the effective temperature on the line profiles in 
subsequent work. However, log(g) = 5.5 as well as [— 1 = 
+0.5 are rather high values and not expected for VB10 (cf. 
e.g. Henry & McCarthy Jr., 1993), so that the lower part of 
the interval is physically more likely. 

Our calculations result in damping constants for alkali 
metals that agree satisfactorily well with the Keck observa- 
tions. The damping constants for non-alkali metals include 
the correction factor of Cg orr = 10 1 ' 8 and lead to line pro- 
files which are consistent with the data. In order to obtain a 
better value for the correction factor, observations of strong 
non-alkali lines are required. We calculated synthetic spec- 
tra outside the observed spectral range available in this pa- 
per in order to find such strong lines, but the only line is 
a moderately strong FelA8829 line besides the RblA7800 
and RblA7948 lines. The iron lines are stronger in hotter 
M dwarfs which would be more suitable for testing the vdW 
broadening of iron lines. We also find that our damping con- 
stants for atomic lines are in broad agreement with the val- 
ues provided by Kurucz. No error should be expected from 
uncertainites in the oscillator stengths of the atomic lines 
since they are much smaller than any other uncertainties 
and are well known for the light element lines that we have 
mostly considered here. 

Our analysis of the Keck spectra yields values for the 
gravity which are consistent with the results of interior cal- 
culations for VB10. The T e g, [— 1 and \og{g) combinations 
agree reasonably well with evolutionary models for M dwarfs 
(Burrows et al., 1993; Baraffe et al., 1995). Larger damping 
constants, e.g. as a result from a different treatment of line 
broadening, would immediately imply a too low log(g). 

There are still discrepancies between models and ob- 
servations in some molecular line dominated regions of the 
spectrum. As mentioned we suspect this to be due to the 
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Figure 8. The fit to the Cal A6573 line. Solid : Observed. Dotted : T eS = 2700K, log(g) = 5.0, [M] = 0.0 (left) and T cff = 2700^, 
log(</) = 5.5, [||] = +0.5 (right). Both models with 8kms _1 rotationally broadened. Note H a in emission. 

Figure 9. The fit to the Bal A7911 and Fcl A7913 lines. Solid : Observed. Dotted : T eS = 2700if, log(p) = 5.0, [M] = 0.0 (left) and 
T cff = 2700A", log(g) = 5.5, [M] = +0.5 (right). Both models with Skins" 1 rotationally broadened. 



molecular line data. This will be improved in subsequent 
work. We also will investigate the influence of Non-LTE ef- 
fects on M dwarf spectra as well as chromospheric features. 
In addition, observations of different strong atomic lines in 
spectral regions where molecular absorption is not dominat- 
ing and blending is less severe (e.g. beyond 1.18/um, see also 
Jones et al., 1996), will clarify the calculations of the damp- 
ing constants and the use of correction factors. 
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Figure 10. The Fel A8662 and Fel A8689 lines. Solid : Observed. Dotted : T cS = 2700^, log(g) = 5.0, [ M ] = 0.0 (left) and T eS = 2700iv , 

log(g) = 5.5, [—J = +0.5 (right). Both models with 8kms _1 rotationally broadened. There has been an offset applied to the observed 
spectrum. See text for details. 



Figure 11. Molecular lines. Solid : Observed. Dotted : T off = 2700i<", log(g) = 5.0, [M] = 0.0 (left) and T eS = 2700K", log(g) = 5.5, 
[— I = +0.5 (right). Both models with 8kms _1 rotationally broadened. 



